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SECTION  I 


INTRODUCTION 

The  shape  of  the  energetic  plasma  spectrum  produced 
when  the  expanding  debris  from  a  high  altitude  nuclear 
detonation  interacts  with  the  background  plasma  has  been 
under  study  and  discussion  for  years.  The  component  of 
the  energetic  particle  distribution ‘parallel  to  the 
magnetic  field  is  particularly  important  to  questions 
concerning  the  location  of  the  ionized  deposition  region, 
the  level  of  ionization  and  the  duration  of  the  ioniza¬ 
tion.  The  answers  to  these  questions  are  necessary  to  do 
an  adequate  job  of  modeling  a  high  altitude  nuclear  event 
(H.A.N.E.)  and  the  subsequent  systems  evaluations  per¬ 
formed  with  these  models. 

Researchers  have  addressed  the  problem  of  field 
aligned  acceleration  of  plasma  in  the  past.  Clark  and 

Papadopoulos1  looked  into  this  problem  in  1976  and 

2 

decided  that  the  Post-Rosenbluth  ,  loss-cone  instability 

would  preferentially  scatter  particles  with  perpendicular 

velocities  of  about  2.5  times  the  thermal  velocity  into 

the  parallel  direction.  This  would  create  a  very 

energetic  spectrum  coming  down  the  field  lines.  They 

addressed  the  nonlinear  three-wave  coupling  of  waves 
.  2 

excited  by  the  Post-Rosenbluth  instability  and  examined 

the  creation  of  a  superthermal  tail  on  the  distribution. 

Then  an  estimate  was  made  of  the  amount  of  energetic 

plasma  scattered  from  the  perpendicular  direction  into 

the  parallel  direction.  Most  of  the  nonlinear  analysis 

was  performed  assuming  that  the  plasma  was  low  beta, 

2 

i.e.,  8  =  8irnkT/B  <1.  A  model  of  this  process  was 

created  and  used  in  the  NRL  codes  for  early  time 


simulations.  It  should  be  noted  that  in  coupling  regions 
of  expanding  debris,  the  plasma  beta  is  expected  to  be 
greater  than  one. 

3 

It  was  later  estimated  by  Palmadesso  that  the 
simple  inverse  mirroring  effects  would  also  help  turn 
particles  from  a  perpendicular  direction  to  a  parallel 
direction  and  would  preferentially  turn  particles  with 
velocities  above  the  debris  streaming  velocity.  This 
work  assumed  that  the  magnetic  moment  of  the  particle 
is  conserved  during  the  transit  down  the  magnetic  field 
line.  This  assumption  could  be  severely  weakened  by  the 
turbulence  expected  to  exist  as  well  as  electromagnetic 
waves  that  might  be  anticipated  in  a  high  beta  system. 

Mission  Research  Corporation  has  also  performed 

4  5 

early  time  H.A.N.E.  calculations  in  the  MHD  limit.  ' 

The  results  from  their  codes  predict  a  much  less 
energetic  spectrum  coming  down  the  field  line.  The 
differences  in  the  NRL  and  MRC  spec t rums  produce 
noticeably  different  results  in  codes  which  do  chemistry 
based  on  the  level  of  ionization  produced  by  plasma 
streaming  down  the  field  lines  and  depositing  in  the 
atmosphere . 

There  are  a  host  of  reasons  for  the  different 
results  and  different  approaches  taken  by  MRC  and  NRL. 
Most  of  them  are  based  on  how  the  respective  groups 
perceived  the  evolution  of  the  interaction  of  the 
streaming  energetic  debris  with  the  background  ambient 
plasma.  The  energy  spectrum  of  plasma  coming  down  the 
magnetic  field  is  a  diagnostic  for  the  larger  picture 
of  this  interaction. 
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If  one  assumes,  as  MRC  does  and  must  because  of  the  MHD 
approach,  that  the  pickup  of  plasma  and  magnetic  field  is  caused 
by  Lamor  coupling,  then  one  is  constrained  to  have  large  scale 
effects.  Any  "model"  for  turning  plasma  transport  and  energy 
into  the  magnetic  field  direction  which  is  based  on  a  fluid 
approach  especially  a  one  fluid  approach  will  produce  a  plasma 
distribution  with  a  maxwellian  shape  centered  about  the  local 
zero  parallel  velocity.  This  occurs  because  of  momentum  and 
energy  conservation  requirements  and  the  iimited  number  of 
degrees  of  functions  in  a  one  fluid  approach. 

The  opposite  view  is  one  where  the  plasma  is  not 
a  fluid  but  exhibits  a  collisionless  behavior.  This 
necessitates  a  different  approach  to  the  dynamics.  Here 
one  might  expect,  wave-particle  and  wave-wave  inter¬ 
actions  with  large  magnetic  field  compressions  occurring 
in  a  much  smaller  interaction  length**'7  than  from  a 
fluid  approach.  The  enhanced  field  compression  and 
small  interaction  region  may  allow  the  magnetic  field 
to  play  a  major  role  in  momentum  conservation  of  the 
particles  being  scattered.  If  the  field  does  take  up 
appreciable  momentum  in  the  scattering  process  that  is 
supposed  to  send  particles  down  the  magnetic  field  lines, 
the  spectrum  of  such  particles  may  hold  little  resem¬ 
blance  to  that  expected  from  the  usual  MHD  approach. 

It  is  implicitly  assumed  by  NRL  in  their  previous  work 
that  this  takes  place.  The  upshot  of  this  preferential 
scattering  is  an  energetic  plasma  system  and  the 
maintenance  of  the  loss  cone  distribution  to  do  the 
scattering.  The  difficulty  here  is  that  in  order  to 
maintain  the  loss  cone  distribution  and  not  end  up  with 
a  lot  of  plasma  moving  in  v^  space  with  the  velocity 
of  the  magnetic  field  lines  one  must  have  a  scattering 
mechanism  that  behaves  as  if  the  scattering  is  reasonably 


7 


inelastic.  This  may  occur  but  probably  only  if  the 

compressed  magnetic  field  region  is  actively  involved 

in  the  scattering  process  such  as  through  electromagnetic 

6  7 

waves.  There  has  been  some  work  by  Clark  et  al.  '  which 
indicates  magnetic  compression  of  the  sort  mentioned 
might  exist.  However,  all  of  this  work  has  been  in  one 
dimension  and  the  1-D  aspect  limits  the  mechanisms  for 
relaxation  of  the  field  compression  as  well  as  explaining 
little  about  the  parallel  dynamics. 

Because  of  this  divergence  in  views  and  the 
importance  of  the  plasma  distribution  coming  down  the 
field  line  to  calculations  of  ionization,  structure/ 
recombination  chemistry  and  emissions  such  as  IR 
emission,  we  have  undertaken  an  investigation  of  whether 
or  not  there  will  be  particle  acceleration  down  the 
field  lines. 

This  report  will  be  limited  to  research  performed 
on  the  loss  cone  type  of  mechanism.  The  research  is 
being  carried  out  with  one  and  two  dimensional  particle 
codes.  The  results  to  be  discussed  in  this  report  are 
from  the  electrostatic  versions  of  these  codes.  We 
will  discuss  not  only  the  results  of  the  simulation 
work  but  the  validity  of  the  model.  In  addition  we  will 
comment  on  the  role  of  the  loss  cone  mode  in  the  H.A.N.E. 
environment  as  well  as  where  it  is  likely  to  be  found. 

The  report  is  organized  into  five  sections. 

Section  II  will  present  our  numerical  methods  and  the 
parameters  used  in  simulations.  In  Section  III  a 
synopsis  of  the  relevant  theory  dealing  with  the  loss 
cone  mode  is  presented.  The  fourth  section  is  devoted 
to  the  results  of  the  simulations  themselves.  In  the 
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last  section  we  will  summarize  and  discuss  our  conclusions 
from  this  research  and  provide  some  guidelines  to  where  it 
should  proceed  from  this  point. 


SECTION  II 


NUMERICAL  METHODS 

The  basic  tool  used  in  this  research  is  the 
particle  code  in  one  and  two  dimensions.  It  differs 
from,  for  instance,  an  MHD  code  in  that  only  the  Lorentz 
force  is  applied  to  the  individual  particles  which 
populate  the  grid.  The  particle  response  is  collision¬ 
less  in  the  classical  sense.  More  importantly,  the 
particle  simulation  allows  kinetic  phenomenon  such  as 
wave-particle  and  wave-wave  interactions  to  occur. 

This  permits  simulation  of  small  scale  (less  than  a 
gyroradius)  collisionless  interaction,  something  that 
by  the  very  assumptions  used  to  obtain  the  MHD  equations 
is  proscribed  in  an  MHD  code. 

The  H.A.N.E.  problem  is  obviously  three 
dimensional  but  because  of  the  usual  computer  limitations 
(size  and  speed)  the  best  we  can  hope  to  achieve  is  a 
2  1/2  dimensional  simulation  (three  velocities  and  two 
spatial  dimensions) .  We  are  interested  in  the  field 
aligned  acceleration  aspect  of  the  total  problem, 
therefore  we  have  selected  the  x  direction  to  be  the 
expansion  direction  across  the  magnetic  field  while  the 
z  direction  is  taken  to  be  along  the  magnetic  field. 

The  usual  loss  cone  analysis  assumes  k( (  <<  k^ 
and  only  addresses  the  quasilinear  effects  of  the 
perpendicular  electric  field.  We  are  looking  for 
acceleration  down  the  field  line  and  must  consider  the 
possibility  that  the  modes  of  interest  may  have  a 
significant  parallel  electric  field.  For  such  fields 
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to  occur  the  parallel  wavelength  shall  be  sufficiently 
long  so  that  electron  Landau  damping  is  negligible 

(to  >>  k  v  )  . 

II  II 

The  requirement  that  X^  <<  X()  produces  a  con¬ 
siderable  hardship  for  a  particle  simulation  because  in 
the  simplest  of  particle  codes  one  cannot  have  cell 
sizes  much  larger  than  the  electron  debye  length  (XD  ) . 
But  to  avoid  the  electron  Landau  damping  one  needs 
waves  in  the  system  from  (10  to  200)  XDe  in  the  parallel 
direction.  This  is  impossible  if  one  is  to  resolve  the 
perpendicular  direction  adequately  enough  to  see  the 
predicted  nonlinear  behavior. 

We  have  greatly  expanded  the  limit  on  the  grid 
size  by  going  to  a  more  complex  interpolation  scheme  in 
the  code,  doing  filtering  of  the  k-space  electric  fields 
and  using  particles  of  macroscopic  size  (larger  than  a 
perpendicular  cell  size) .  The  interpolation  scheme 
selected  and  the  particle  shape  used  is  that  of  a  cubic 
spline.  The  actual  shape  function  has  the  following 
form: 


S3  (x) 


x2(2  -  I  x  | )/  2^1 


2/3 

(2  -  j  x |  )  3/6 


.1  < 


M  i  i 
M  <  2 

|x|  >  2 


(1) 


where  x  denotes  the  displacement  of  the  individual 
particles  from  the  surrounding  grid  points.  Basically 
the  cubic  spline  interpolation  is  used  to  distribute 
the  particle  charge  over  the  grid  and  to  apportion  the 
resultant  electric  fields  to  the  particle  positions.  A 
cubic  spline  Fourier-space  particle  shape  factor  is 
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applied  to  smoothly  cut  off  the  k-spectrum  and  preserve 
integrity  in  real  space.  A  smooth  cut-off  of  the 
spectrum  is  important  since  much  of  the  numerical  grid 
noise  is  carried  by  wavelengths  on  the  order  of  the 
grid  spacing  (largest  k-modes) .  An  abrupt  cut-off  with 

fl  I  u  1  <2 


u  |  >  2 


where  u  =  k/(l/2  k  ) ,  would  yield  a  real-space 

IHclX 

particle  shape  in  the  form  of  a  sine  function 


iin(1/2  rW/(1/2  rkmax> 


which  has  the  undesirable 


property  of  negative  "ringing"  attenuated  only  as  1/r. 

4  4 

Cubic  splines  go  as  sine  and  produce  a  1/r  attenuation. 

The  result  of  the  cubic  spline  approach  which  is  costly 

on  a  per  particle  basis  is  that  we  can  use  grid  cells 

of  10  to  40  XDe  in  the  parallel  direction.  We  have  in 

fact  done  simulation  tests  to  make  sure  that  the  grid 

effect  problem  has  been  reduced  to  a  negligible  level, 

and  have  found  this  to  be  the  case. 

The  one  and  two  and  a  half  dimension  simulations 

reported  here  are  electrostatic  in  nature.  In  the  1-D 

code  we  used  a  256  cell  mesh  where  each  cell  was  2  A_  . 

De 

A  series  of  tests  were  run  looking  at  the  behavior  of  a 
counter  streaming  ion  beams  using  1,  4,  8,  16  particles  per 
cell.  The  purpose  was  to  determine  how  many  particles/ 
cell  were  necessary  for  an  adequate  simulation.  The 
effect  of  too  few  particles  can  be  seen  in  Figures  (1-2) 
where  in  addition  to  the  relaxation  of  the  distribution  to 
a  more  maxwellian  shape  particle  acceleration  was  also 
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Figure  1.  A  histogram  of  the  ion  distribution  function 

for  a  run  of  4  particles/cell.  Note  scattering 
to  high  positive  and  negative  velocities,  this 
occurs  because  of  numerical  noise  in  the  system. 
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Figure  2.  A  histogram  of  the  ion  energy  flux  for  a  run  of 
4  particles/cell.  This  shows  the  large  change 
in  energy  flux  caused  by  the  numerically 
generated  diffusion  in  the  simulation. 


seen  in  the  high  energy  tail  of  the  distribution  due  to 

spurious  electric  fields.  At  8  particles/cell  the 

numerical  effects  are  greatly  reduced  (Figures  3-4)  and 

at  16  particles/cell  (Figures  5-6),  the  whole  simulation 

was  extremely  quiet.  For  the  work  to  be  discussed  in 

this  report  8  particles/cell  were  selected. 

The  2  1/2-D  code  used  an  x,  z  grid  of  128  by  32 

cells  where  in  the  x  direction  the  cell  size  is  2  A_ 

De 

and  in  the  parallel  direction  the  cell  size  is  40  A  . 

In  both  codes  the  nominal  mass  ratio  m./m  is  64/1. 

i  e 

For  the  actual  H.A.N.E.  simulation  it  was  assumed  that 
the  loss  cone  distribution  would  be  maintained  as 
discussed  in  Clark  and  Papadopoulos. 1  Therefore  an 
algorithm  was  placed  in  the  code  which  tends  to  maintain 
the  loss  cone.  The  importance  of  maintaining  this  form 
will  be  discussed  in  the  next  section.  Also  a  review 
of  theory  will  be  presented  in  that  section. 
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Figure  3. 


A  histogram  of  the  ion  distribution  function 
for  a  simulation  with  8  particles/cell.  The 
noise  has  been  substantially  reduced  and  much 
less  scattering  occurs  particularly  in  the 
larger  maxwellian  distribution. 
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Figure  4.  A  histogram  of  the  ion  energy  flux  for  8  particles/ 
cell.  The  amount  of  diffusion  to  high  energies 
has  been  significantly  reduced  but  is  still  in 
evidence . 
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Figure  5.  A  histogram  of  the  ion  distribution  function 
for  a  simulation  with  16  particles/cell.  It 
is  quieter  than  is  seen  in  Figure  4  but  only 
marginally. 
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Figure  6.  The  histogram  of  the  ion  energy  flux  for  16 

particles/cell.  The  amount  of  diffusion  seen 
here  is  essentially  the  same  as  Figure  4. 


SECTION  III 


THEORY 

The  term  loss  cone  instability  is  actually 

applied  to  a  larger  class  of  instability  than  is  found 

in  a  magnetic  mirror  machine.  In  the  DNA  community  the 

term  loss  cone  also  refers  to  a  physical  region  in  the 

perturbed  magnetic  field  where  plasma  simply  streams 

down  the  magnetic  field  lines.  The  loss  cone  instability 

is  often  applied  to  situations  where  there  is  a  deficiency 

in  the  parallel  part  of  the  distribution  through  loss  or 

enhancement  of  a  portion  of  the  perpendicular  velocity 

space  such  that  a  gradient  in  the  distribution  f (v^ ,  v() 

is  produced  (Figure  7) .  Unstable  loss  cone  distributions 
„f  ion 

- (v  ,  v  )  >0,  are  found  in  the  fusion  devices  such 

J_  ^  ^  g  g 

as  magnetic  mirror  machines  (Gerver  and  Cohen  )  and 
tokomaks  where  neutral  beam  injection  is  taking  place, 
Hitchcock  et  al.10  They  exist  in  the  ionosphere  particu¬ 
larly  in  the  polar  region  and  in  the  magnetosphere.  This 
instability  is  predicted  to  play  a  critical  role  in  the 

Southern  Conjugate  rebrightening  seen  in  Starfish, 

11  12 

Brecht  et  al.  '  .  Of  most  importance  to  this  work, 

the  loss  cone  instability  was  predicted  by  Clark  and 
Papadopoulos1  to  produce  in  its  nonlinear  evolution  an 
acceleration,  or  perhaps  scattering  is  a  better  word,  of 
high  energy  plasma  down  the  magnetic  field  line  during 
a  H.A.N.E. 

The  linear  analysis  that  is  traditionally  per¬ 
formed  assumes  that  L,  the  macroscopic  scale  length, 
i.e.,  the  length  of  the  field  line,  is  large  compared  to 
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Loss  cone - — 


Figure  7. 


a)  This  is  the  typical  distribution  found  in 
mirror  geometries.  This  distribution 
occurs  in  mirror  machines  and  at  the  ends 
of  geomagnetic  field  lines. 

b)  This  is  a  distribution  often  termed  loss  cone 
produced  by  an  unmagnetized  plasma  beam 
interacting  with  a  maxwellian  background. 

c)  This  is  the  loss  cone  region  often  referred 
to  in  H.A.N.E.  work.  This  geometry  is  the 
same  as  the  mirror  machines. 


the  wavelength,  A  ,  of  the  loss  cone.  This  assuption, 

L  >  A  ,  is  made  so  that  the  mode  has  sufficient  room 
II 

to  grow  as  it  propagates  down  the  magnetic  field 

line.  A  second  assumption  is  that  >>  k  (  where 

k  =  2n/A.  Also  it  is  generally  assumed  that  the 

frequency  of  the  wave,  u  ,  is  much  greater  than  the 

^  eB 

ion  cyclotron  frequency,  co  .  =  —  and  that 

2  2  me 

k^  pf  >>  1  where  is  the  ion  gyroradius.  This 

assumption  means  that  to  first  order  at  least  the  ions 

are  unmagnetized;  a  good  assumption  in  the  early  time 

H.A.N.E. 


A  further  assumption  about  the  electrons  is 

2  2 


usually  made  and  that  is  |  |  <<: 


(the  electrons  are  magnetized) . 
relaxed  and  will  need  to  be  in  the  H.A.N.E 


a)  and  k  p  <  <  1 

1  ce 1  1  e 

This  assumption  can  be 


case . 


The 


relaxation  of  this  assumption  only  reduces  the  growth 

2  2  1/2 

rate  by  1/(1  +  w  /w  )  '  where  w  is  the  plasma 


pe  ce 
.2.1/2 


pe 


frequency,  =  •  In  its  simplest  form  the  real 

frequency  of  the  eioss  cone  instability  for  high 
velocities  is  found  to  be 


to 


k 


(3) 


with  a  growth  rate  that  is  approximately 


(4) 


where  6(w/kv)  is  a  sinusoidal  function  with  argument 

(w/kv  -  it)  and  v  is  the  mean  thermal  speed. 

13  14 

As  shown  by  Davidson  and  Galeev  ,  the  nonlinear 
evolution  of  the  loss  cone  instability  diffuses  particles 
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into  the  loss  cone  region  and  thus  stabilizes  the  wave  by 
3F 

reducing  ■*  0.  However,  if  the  loss  time,  t^,  out 

of  the  system,  i.e.,  down  the  field  line,  or  plasma 
replenishment  time,  tr»  is  such  that  cold  plasma  in  the 
moving  magnetic  field  frame  of  reference  in  the  H.A.N.E. 
case  is  interacting  with  the  system  faster  than  the  non¬ 
linear  diffusion  time,  tD»  the  loss  cone  instability  will 

persist;  t  or  t  <  t_.  If  such  a  situation  exists  then, 
r  x,  D  15 

as  pointed  out  by  Galeev  ,  the  only  way  for  the  unstable 

modes  to  limit  themselves  is  through  mode  coupling  to  a  mode 

that  has  a  plasma  velocity,  w/k,  that  matches  a  portion  of 

3F 

the  distribution  —  <  0  where  it  will  damp  and  place 

the  energy  back  into  the  distribution.  This  produces  a 
diffusion  of  the  negative  slope  of  the  ion  distributions 
to  high  energy.  Clark  and  Papadoupoulos1  pointed  out 
that  such  a  case  could  exist  for  a  H.A.N.E.  They 
theorized  that  if  this  nonlinear  state  were  reached  then 
particles  with  high  velocity  would  be  scattered  off  of 
the  high  phase  velocity  waves  down  the  magnetic  field 
line.  They  arrived  at  an  estimate  of  2.5  vfchi  as  the  average 
scattering  velocity,  where  vfc^  is  the  ion  thermal  velocity. 

However,  if  one  considers  the  type  of  plasma 
parameters  predicted  by  the  multifluid  code,  KLYSMA6, 
in  the  interaction  region,  one  finds  that  for  oxygen 
or  aluminum  that  vthi  at  15  keV  is  approximately 
2.5-5  times  less  than  the  expected  2  x  io®  streaming 
velocities  and  therefore  the  modes  driven  unstable  by 
the  loss  cone  could  provide  the  scattering  parallel  to 
the  magnetic  field  line  but  would  not  provide  particles 
moving  faster  than  the  original  expansion  velocity. 

However,  the  estimated  phase  velocity  of  these  waves 


parallel  to  the  field  line  does  approach  the  electron 
thermal  velocity.  Therefore,  scattering  by  such  waves 
may  produce  velocities  approaching  the  parallel  phase 
velocity  and  be  even  more  energetic  than  predicted  by 
Clark  and  Papadopoulos . * 

One  of  the  interesting  aspects  of  the  current 
loss  cone  mode  analysis  is  that  most  of  the  nonlinear 
treatments  only  consider  the  perpendicular  diffusion 
coefficient.  However,  the  existence  of  a  finite  k|( 
is  in  fact  postulated  and  used  to  estimate  growth  along 
the  magnetic  field  line.  Therefore  one  in  fact  has  an 
E^  associated  with  the  mode.  This  E^  driven  to 
sufficiently  high  levels  by  nonlinear  waves  could 
provide  the  predicted  scattering. 

All  of  this  discussion  leads  to  the  basic  reason 
we  are  performing  this  research.  There  appears  to  be 
several  possible  mechanisms  by  which  the  loss  cone 
instability  might  provide  an  energetic  field  aligned 
spectrum.  The  pure  scattering  mentioned  by  Clark  and 
Papadopoulos,  nonlinear  energy  diffusion  in  the  mode 
spectrum  allowing  parallel  propagation  waves  to  scatter 
particles,  or  electromagnetic  responses  due  to  inhomo¬ 
geneous  magnetic  fields  are  all  candidates.  If  any  of 
these  mechanisms  work,  it  is  through  a  very  nonlinear 
interaction  between  various  waves  and  the  distribution 
function.  The  analysis  drawn  from  data  is  only  an 
approximation,  therefore  in  order  to  determine  what 
actually  does  occur  when  a  loss  cone  distribution  is 
driven  to  a  nonlinear  state  we  have  chosen  to  simulate 
it  with  a  particle  code. 

In  the  next  section  we  will  discuss  these  simula¬ 
tions  performed  to  date  and  their  results. 


SECTION  IV 


SIMULATION 

During  the  course  of  the  contract  two  sets  of 
simulation  codes  were  utilized;  a  1-D  cubic  spline  code 
and  a  2  1/2-D  cubic  spline  code.  The  1-D  code  was  built 
and  used  as  a  test  bed  to  develop  numerical  algorithms 
and  to  determine  appropriate  physical  parameters.  Once 
the  appropriate  determinations  were  made  for  numerical 
and  physical  parameters  the  2  1/2-D  code  was  used  to 
study  the  loss  cone  problem. 

The  1-D  code  was  run  with  256  cells  and  8  particles/ 
cell.  The  cell  size  decided  upon  was  Ax  =  2  ^De»  All 
of  the  tests  were  run  with  the  one  dimension  perpendicular 
to  the  magnetic  field.  The  initial  distributions  for  these 
runs  were  maxwellian  electrons  and  ion  distributions  of  the 
form  shown  in  Figures  1,  3  or  5.  Initially  fusion  type 
parameters  were  run  with  both  ions  and  electrons  magnetized. 
Later  only  the  electrons  were  magnetized,  ajce/Cl)pe  %  2r 
and  the  ions  free  streamed  consistent  with  the  H.A.N.E. 
situation.  The  ions  were  loaded  onto  the  simulation  grid 
as  two  maxwellians  of  10  keV  temperature.  One  quarter 
were  loaded  into  one  distribution  and  three  quarters  into 
the  other.  A  relative  drift  of  five  times  the  ion  thermal 
velocity  was  imparted  to  the  test  distribution  similar  to 
possible  H.A.N.E.  situations.  This  gave  a  center  of  mass 
velocity  of  3.75  times  the  thermal  velocity.  As  can  be 
seen  in  Figure  3  the  initial  distribution  was  not  that  of 
a  classic  loss  cone  in  v^  space,  but  was  more  the  beam 
plasma  type  of  distribution.  After  t  ^  10/wpe  the 
diffusion  caused  by  the  instabilities  in  the  plasma  have 
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created  a  loss  cone  like  distribution  as  seen  from  the 
response  of  non-drifting  distribution.  From  here  on 
the  evolution  should  follow  that  of  a  loss  cone  distribu¬ 
tion.  One  important  point  to  be  made  here  with  regard 
to  all  these  distributions  is  that  the  classic  loss  cone 
analysis  assumes  a  distribution  function  which  is 
symmetric  about  the  magnetic  field.  Until  one  reaches 
times  of  'vO.Ol  second  or  longer  after  the  debris  and 
air  interaction  begins  such  assumptions  for  the  H.A.N.E. 
case  (Starfish)  are  not  appropriate.  More  discussion  on 
this  topic  is  forthcoming  in  the  conclusion  section  of 
this  report. 

A  variety  of  mundane  but  important  parameters 
were  tested  with  the  1-D  code.  The  effect  of  the  ion/ 
electron  mass  ratio  was  tested  by  varying  the  ratio  from 
16/1  to  128/1  in  powers  of  two.  This  showed,  as  expected 
from  theory,  that  the  effects  seen  in  the  simulation  were 
caused  by  the  ions  primarily  although  the  electrons  do 
play  a  role  in  the  dynamics .  The  mass  ratio  decided 
upon  for  the  major  simulation  was,  64/1,  ion  to  electron 
ratio.  The  effect  of  grid  size  was  also  checked  to 
determine  if  the  size,  Ax,  of  the  grid  provided  both 
sufficient  spectral  density  and  range  of  k  spectrum  to 
allow  the  expected  dynamics  such  as  resonant  interaction 
with  the  slope  of  the  distribution  and  three  wave  mode 


coupling.  The  grid  of  Ax  ^  2  ADe  was  found  to  be 
sufficient  to  represent  the  appropriate  resonance 
regions  as  well  as  allow  the  nonlinear  mode  coupling 


expected;  ADe  is  the  electron  debye  length 

ADe  =  vtjle/(i0pe*  An  important  result  was  that  one  could 

use  a  system  of  128  cells  in  the  perpendicular  direction 
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and  perhaps  64  if  necessary  and  get  adequate  resolution 
with  a  periodic  system.  This  information  was  particu¬ 
larly  important  when  the  2  1/2-D  code  was  to  be  used 
because  we  then  knew  that  128  cells  in  a  doubly 
periodic  system  would  be  sufficient  leaving  some  memory 
and  money  to  include  the  second  dimension  and  required 
number  of  particles. 

Figure  3  demonstrates  the  dynamic  effect  the 

plasma  instabilities  have  on  the  distribution  function. 

The  total  time  covered  in  Figure  3  is  160/uj^  where 

each  line  drawn  represents  At  =  5/to  .  As  can  be  seen 

pe 

from  Figure  3  diffusion  to  high  velocities  occurred 
for  both  distributions,  however,  the  lower  density 
distribution  is  seen  to  be  completely  smoothed  out  by 
the  end  of  the  run  and  bears  no  resemblance  to  a 
maxwellian.  This  occurs  because  of  the  center  of  mass 
effect.  This  plot  is  shown  in  the  systems  center  of 
mass  which  explains  why  the  larger  density  distribution 
function  is  slightly  diffused.  Note  that  by  t  ^  20/uj^ 
strong  velocity  diffusion  has  almost  filled  in  the 
hole  between  the  two  distributions.  Figure  4  shows  the 
effect  of  this  diffusion  on  the  energy  flux  of  the 
system.  Here  one  sees  a  rather  substantial  shift  upward  in 
the  energy  flux.  It  should  again  be  emphasized  that 
this  was  a  very  energetic  beam  distribution  with  two 
distinct  ion  species  not  a  thermalized  loss-cone  dis¬ 
tribution  with  one  ion  species. 

The  electric  field  spectrum  behaved  as  expected. 

The  total  electric  field  energy  was  approximately  ten 
percent  of  the  thermal  energy  in  the  system.  The  mode, 
k  =  25,  that  resonated  with  the  positive  slope  of  the 
small  beam  was  excited  to  high  levels  as  compared  to 


all  other  modes.  This  can  be  seen  in  Figure  8  which 
shows  all  the  modes  in  the  system  including  the  con¬ 
jugate  ones.  Because  of  the  periodic  boundary  conditions 
the  positive  and  negative  modes  are  identical.  In  this 
section  the  dimensional  mode,  K.,  is  defined  as 

_  •  i  1 

K.  =  -r— rr  =  -r—  k ,  where  N  is  the  number  of  cells  in  a 
specific  direction  and  k  is  the  dimensionless  mode 
number  referred  to  in  this  section  which  corresponds  to  the 
index,  i,  i.e.,  k  =  25.  Also  seen  in  Figure  8  is  that 
k  =  2,  a  long  wavelength  mode  with  a  phase  velocity 
that  matches  particle  velocities  far  out  in  the  distri¬ 
bution,  is  the  only  other  dominant  mode.  These  two  modes 
vary  dramatically  during  the  simulation  as  can  be  seen 
in  Figures  9  and  10.  However,  only  the  k  =  25  mode  is 
obtaining  its  energy  directly  from  the  beam.  The  long 
wavelength  mode  is  apparently  deriving  its  energy  from 
nonlinear  three  wave  coupling. 

Looking  at  Figure  10,  one  can  estimate  the  linear 
growth  rate  of  mode  25  assuming  an  exponential  growth. 

Making  this  assumption  a  growth  rate  of  y  ^  . 12  w  is 

16  pe 

obtained.  Using  a  maximum  entropy  technique  the 

fundamental  frequency  of  this  mode  was  determined  to  be 

^  .08  u  .  The  rate  of  diffusion  and  growth  for  the 
pe 

modes  can  be  checked  with  theory  to  determine  if  the 

simulation  is  behaving  as  quasilinear  theory  might 

suggest.  Although  most  of  the  theory  is  for  a  two 

dimensional  system  the  quasilinear  diffusion  is  generally 

a  one  dimensional  calculation.  Using  the  work  of 
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Sagdeev  and  Galeev  we  employ  their  equation  for  quasi¬ 
linear  diffusion  (11-59) 
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Figure  9.  Time  evolution  of  mode  two,  k 
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Figure  10.  Time  evolution  of  mode  twenty  five 
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where  w  =  v  /vtj^  #  y  =  vt^)  and  i|>  is  the 

distribution  function.  One  can  estimate  the  time 
necessary  for  quasilinear  diffusion,  t^,  from  Eq.  (5) 
by  making  the  following  assumptions: 


ii  %  Jl 
at  t 

q 


and 


%  “k/ki 


using  these  assumptions  the  following  relation  is 
obtained: 


ve 


i*l 


~~2  3  T~.  ..1/2  ..  , 

m^  (Aw) 
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A  reasonable  assumption  for  Av,  is  that 


„  2 


thi 


With  this  assumption,  the  information  from  Figure  10 
and  the  real  frequency,  x^  can  be  evaluated,  where  we 
have  correctly  assumed  that  k  =  25  produces  most  of 
the  diffusion  effects.  Using  this  information  one 
obtains  the  result  that 


'v  (.36  u)  . ) 
q  pi 


-1 


This  implies  that  the  quasilinear  diffusion  time  varies 
1/2 

as  nr  '  .  This  was  checked  against  the  runs  made  with 
varying  mass  ratio  and  found  to  agree  quite  well.  In 
terms  of  the  mode  data  from  the  8  particles/cell 
simulation,  x^  ^  2®/wpe  usi-n9  the  64/1  mass  ratio. 
This  is  in  very  good  agreement  with  the  data.  Figure  3. 

In  order  to  determine  if  nonlinear  three  wave 
interaction  took  place  further  work  with  the  electric 


fields  was  necessary.  Figures  9-11  show  the  time  history 
of  the  most  important  modes.  Figure  9  shows  the  time 
history  of  mode  two.  Note  the  rapid  growth  in  the  energy 
of  the  mode.  Figure  10  displays  the  main  resonant  mode, 
k  =  25.  As  inferred  from  Figure  8,  this  mode  has  most 
of  the  electric  field  energy.  Figure  11  shows  mode  27. 

The  data  shown  on  these  plots,  as  well  as  the  other 
modes,  were  analyzed  to  determine  their  frequencies.  A 
standard  FFT  was  found  to  be  inadequate  for  this  task 
because  of  the  sparsity  of  temporal  data.  Instead  we 
employed  a  maximum  entropy  technique"^  which  is  consider¬ 
ably  more  sensitive  to  sinusodial  variation  in  a  sequence 
of  data.  It  was  presumed  in  this  analysis  that  a 
sinusoidal  dependence  existed  and  certainly  Figures  9-11 
demonstrated  that  it  does  for  the  most  part.  Because 
of  this  assumption  and  the  fact  that  one  uses  coefficients 
to  determine  the  frequency  curve, spurious  signals  can  be 
generated.  We  used  both  a  5  and  a  25  coefficient  curve 
to  check  our  results.  There  were  differences  in  the 
results.  Only  the  frequency  with  the  largest  PSD  for 
each  mode  was  considered.  One  could  reasonably  expect 
to  see  several  other  modes  besides  the  driven  ones  in 
the  PSD's  and  they  do  appear.  The  electron  plasma 
oscillation  should  be  present, as  well  as, the  ion  plasma 


frequency,  .25  w  ,  and  the  electron  sound  wave  wil 

2  2  2  Pe 

a)  %  w  +  to  Evidence  for  all  these  modes  was  foi 
pe  ce 

Using  the  maximum  entropy  technique  we  obtained 
the  required  data  to  check  and  see  if  nonlinear  wave 


and  the  electron  sound  wave  with 


Evidence  for  all  these  modes  was  found. 


interactions  were  occurring.  Specifically  we  looked 
for  waves  that  met  the  following  criterion  for  wave-wave 
interaction: 
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Figure  11.  Time  evolution  of  mode  twenty  seven,  k  =  27, 
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where  k  and  oj  may  be  either  positive  or  negative. 
Looking  at  a  series  of  plots  like  those  in  Figure  8 
gave  evidence  that  mode  2  and  25  might  be  coupled 
energetically.  In  Figures  12,  13  and  14  the  frequency 
dependence  of  modes  2,  25  and  27  are  displayed.  It  was 
determined  that  these  modes,  2,  25,  and  27,  met  the 
appropriate  conditions  for  a  three  wave  interaction, 
k2^  +  k_25  =  ^2  and  =  co2.  More  important 

27  and  25  grew  fast  enough  to  account  for  the  rapid 
onset  of  mode  2  as  seen  in  Figure  9.  Mode  27  was  also 
found  to  be  driven  by  the  beam  distribution  but  trans¬ 
ferred  its  energy  into  mode  2.  Other  three  and  four 
wave  interactions  probably  occur  in  this  simulation  but 
because  of  the  rather  monochromatic  behavior  of  modes  2 
and  25  other  relations  were  not  considered  here.  The 
main  purpose  was  after  all  to  determine  if  resonant 
and  nonlinear  behavior  could  occur  in  the  spectral 
spacing  (Ax)  and  spread  (NAx)  given.  Looking  at 
Figure  8  one  also  notes  that  most  of  the  dynamics  in 
the  electric  field  occurred  in  modes  2-35  implying  that 
a  256  cell  system  with  128  modes  and  their  complex 
conjugates  could  be  easily  handled  with  a  128  system 
and  64  unique  modes. 

An  additional  set  of  runs  were  made  where  all 
conditions  discussed  above  were  repeated  except,  that 
instead  of  letting  the  distribution  relax  normally,  a 
source  of  plasma  was  placed  into  the  simulation  to 
represent  the  effect  of  cold  plasma  being  constantly 
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Figure  12.  Power  spectral  density  of  mode,  k.  =  2 


swept  up.  To  do  this  an  algorithm  was  developed  that 
attempted  to  hold  the  initial  distribution  fixed.  Simply 
putting  particles  at  a  specific  velocity  will  lead  to  a 
numerical  beaming  instability.  Instead  a  random  number 
generator  was  used  to  select  a  given  number  of  particles 
from  the  total  velocity  distribution  and  randomly 
redistribute  them  in  velocity  space  so  as  to  match  the 
original  distribution.  For  the  additional  simulations 
approximately  2  particles/step  were  replaced.  This 
represents  a  rather  high  energy  replacement  compared  to 
the  real  situation.  The  results  were  found  to  be  quite 
similar  particularly  in  the  electric  fields.  As  can  be 
seen  by  comparing  Figure  3  with  Figure  15  there  is  a 
slight  change  in  the  late  times  with  the  replenished 
system  maintaining  a  sharper  peak  in  the  beam  and  pro¬ 
ducing  a  sharper  PSD. 

The  2  1/2-D  simulation  was  carried  out  on  a 
standard  loss  cone  distribution.  This  distribution  differs 
from  the  one  dimensional  simulation  as  can  be  seen  in 
Figure  16.  We  studied  the  standard  loss  cone  because  it 
was  specifically  discussed  in  the  literature  as  a 
possible  scattering  mechanism.  The  results  of  the  1-D 
runs  would  imply  a  beam  plasma  interaction  as  a  likely 
condidate.  However,  for  thoroughness  and  completeness 
we  persisted  with  a  thermalized  loss  cone  distribution. 

In  this  case  there  really  is  only  one  ion  distribution 
function.  The  simulation  was  prepared  by  loading  ions 
and  electrons  to  represent  a  5  keV  temperature  and  an 
effective  mirror  ratio  of  two.  This  was  accomplished 
with  a  random  number  generator  as  were  the  1-D  simula¬ 
tions.  The  equation  used  for  loading  was: 
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Figure  15.  Histogram  of  the  ion  distribution  function 

for  a  simulation  where  the  initial  distribution 
is  being  replenished. 
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Figure  16.  The  initial  ion  loss  cone  distribution  for  the 
2  1/2-D  runs.  A  mirror  ratio  of  two  was 
assumed.  This  provided  a  velocity  peak  at 
approximately  2  xlO®  cm/sec. 
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where  R  is  the  mirror  ratio,  a  =  m/2kT,  and  T  is 
the  parallel  or  perpendicular  temperature. 

The  choice  of  temperature  and  mirror  ratio  were 
made  to  be  consistent  with  the  energy  available  to  the 
counterstreaming  ions  in  the  Starfish  case.  Here  we 
attempted  to  emulate  a  loss  cone  distribution  in  per¬ 
pendicular  velocity  space  that  had  a  peak  in  it  at  about 
0 

1  to  2  x  10  cm/sec  where  the  fact  that  the  ions  in  the 
simulation  are  light  was  considered. 

As  discussed  earlier,  it  was  decided  after  testing 
that  using  Az  =  40  XDe  would  be  stable  and  with  appro¬ 
priate  use  of  the  splines  finite  grid  effects  common  to 
such  large  cell  sizes  would  be  overcome.  This  was 
confirmed  in  tests  where  we  detected  no  spurious  heating 
of  the  particles,  particularly  the  electrons,  energy  did 
not  build  up  in  the  shortest  parallel  wavelength  modes 

and  we  conserved  energy  to  the  same  accuracy  as  the  1-D 

-4 

simulations,  (E  -  E  )/E  ^  10 

The  evolution  of  the  loss  cone  distribution  can  be 
seen  in  Figures  16-18  where  the  contours  along  the  v^  =  0 
line  are  seen  to  expand  outward  smoothing  the  contours.  As 
expected  from  quasilinear  theory  the  loss-cone  fills  in  the 
depletion  in  v™  space  but  it  does  not  appear  to  diffuse  to 
high  energy  in  v^  space  as  the  1-D  beam  case  did.  At  the 
time  shown  in  Figure  18  the  process  is  not  complete  but  there 
is  no  evidence  to  suggest  that  nonlinear  diffusion  in  the 
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perpendicular  velocity  space  will  occur  to  any  appreciable 
extent  at  this  late  time  in  the  simulation. 

The  electric  field  behavior  was  very  interesting. 

The  total  electric  field  energy  in  the  system  saturated 
at  approximately  7  *  10  ^  of  the  thermal  energy.  This 
was  a  factor  of  at  least  ten  less  than  the  one  dimensional 
simulation  with  the  same  number  of  particles  per  cell,  8, 
but  no  extra  drift  velocity  energy  available  to  the  system. 
There  was  a  very  definite  k()  dependence  to  the  per¬ 
pendicular  electric  field,  E.  .  The  parallel  electric 

K1 

field,  E,  ,  demonstrated  both  a  k  and  k  dependence 
Kll  1 

but  eventually  the  mode  with  the  most  energy  was  one  where 

A  =  a  •  It  was  also  noted  that  the  peak  parallel 
1  II 

potential,  ,  was  approximately  one  half  of  the 

perpendicular  peak,  for  different  modes.  They  did 

not  peak  at  the  same  time.  Figures  19  and  20  are 
representative  of  the  type  of  dependences  that  were  found. 

The  mode  that  grows  earliest  is  one  where 
A^/A()  ^  1/40  consistent  with  most  analytical  work.  As 
the  simulation  progresses,  the  energy  in  the  perpendicular 
electric  field  cascades  into  the  flute  mode  where  k|(  =  0. 
Ultimately  the  preponderance  of  energy  is  spread  among 
modes  having  one  of  four  available  parallel  components. 
Figures  21,  22,  and  23  show  the  temporal  variations  in 
the  whole  electric  field  spectrum.  There  is  little  doubt 
that  some  nonlinear  wave-wave  coupling  is  occurring  but 
not  to  the  extent  seen  in  the  one-dimensional  simulations. 
As  can  be  seen  by  comparing  e$/kT  for  both  sets  of 
simulations,  Figures  (8  and  19)  the  2  1/2-D  runs  produced 
much  lower  levels  of  potential  fluctuation  per  mode. 
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Figure  17.  The  ion  loss  cone  distribution  at  t  %  80/w 
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Figure  19.  Shown  here  is  the  mode  spectrum  for  the 

perpendicular  electric  field  energy.  All  of 
the  spectrum  is  shown  for  k|(  ranging 

from  one  to  eight.  The  maximum  ei|>/kT  shown 

here  is  9.4  x  lo  Two  sets  of  data  are 

plotted  on  each  ordinate  axis  ranging  from  1  to 
64  and  then  repeated. 
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Figure  21.  Perspective  plots  of  the  electric  field  energy 

at  t  =  9.2/  .  after  the  initial  linear 

pe 

growth  of  the  instabilities.  (a)  The  perpen¬ 
dicular  field  energy.  (b)  The  parallel  field 
energy. 
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Figure  22.  The  electric  field  energy  at  t  =  80.2/oj  . 

Figures  21-23  are  all  self  normalized  p 
comparison  can  only  be  made  by  noting  the 
maximum  amplitude  on  each  plot.  (a)  The 
perpendicular  field  energy.  (b)  The  parallel 
field  energy. 
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The  frequency  for  the  largest  amplitude  mode  in 

the  system,  k  =  (k  ,  k  )=  (4,17),  should  from  equation 

1  111 

(3)  be  0)  *\»  -jpjr-  0)  'V  .03  ai  .  Estimates  from  the  simula- 
40  pe  pe 

tions.  Figure  24,  produce  y  ^  .11  and  u>  'v  .04(0^ 

considering  finite  particle  width  effects  of  dispersion 

relations.  Using  e<j>/kT,  the  wavelength  of  the  modes 
2  2 

and  a  Av^  'u  (vth^/4)  one  can  estimate  the  quasilinear 
diffusion  time  from  Eq.  6.  The  result  is  that  t  ^  56/to  .. 

q  pi 

This  agrees  with  our  result  reasonably  well.  Figures  16 
and  17. 

The  2-D  simulation  was  performed  with  a  replenish¬ 
ment  algorithm  in  it.  Here  we  replenished  at  the  rate  of 
10  particles  randomly  throughout  the  distribution.  One 
can  estimate  what  the  particle  flux  into  a  simulation 
region  such  as  used  here  would  be  in  the  H.A.N.E.  case. 

O 

For  the  estimate  I  chose  v  ^  1  *  10  cm/sec  and 
8  3 

nQ  ^  7  x  io  particle/cm  .  This  resulted  in  a  flux  into 
the  grid  of  0.22  particles/w^  where  for  the  simulation 
At  =  .2/<jjpe.  Because  of  the  attempt  by  the  algorithm  to 
maintain  the  original  distribution  we  feel  that  the  effect 
of  incoming  plasma  to  the  interaction  region  is  more  than 
adequately  modeled  here. 

In  summary  we  have  begun  research  into  possible 
field  aligned  acceleration.  One  dimensional  simulations 
were  performed  using  beams  of  unequal  density.  This  is 
throught  to  be  the  case  for  regions  of  the  interactions 
in  the  H.A.N.E.  It  was  found  that  large  potentials  were 
generated  in  very  discrete  modes.  Strong  diffusion  was 
produced  in  a  fashion  consistent  with  theory  and  earlier 
simulations  for  counterstreaming  plasmas.  The  less  dense 
of  the  two  beams  experienced  energy  diffusion  both  up  and 
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down  the  energy  spectrum  while  the  dense  beam  experienced 
a  small  amount  of  high  energy  diffusion.  The  diffusion 
was  very  fast,  occurring  on  the  time  scale  of  twenty 
electron  plasma  periods.  In  the  two  dimensional  simula¬ 
tions  a  true  loss  cone  simulation  was  undertaken.  The 
electric  field  energy  was  found  to  be  considerably  lower 
than  for  the  1-D  simulation.  In  addition  the  energy  in 
the  perpendicular  electric  field  was  shared  by  many  modes 

and  had  a  k  dependence.  A  parallel  electric  field  was 

1 1 

produced.  The  maximum  potential  found  in  a  parallel  mode 
was  approximately  one  half  of  the  maximum  potential  found 
in  a  perpendicular  mode.  The  level  of  maximum  potential 
for  any  given  mode  was  orders  of  magnitude  below  that  found 
in  the  1-D  case.  As  expected  the  loss  cone  filled  but 
little  diffusion  to  high  energy  was  noted  in  the  perpen¬ 
dicular  portion  of  the  distribution. 

In  the  next  section  a  discussion  of  the  implications 
of  these  results  and  the  questions  they  raise  will  be 
presented.  In  addition  ideas  about  future  directions  will 
also  be  presented. 


SECTION  V 


CONCLUSIONS 

The  results  discussed  in  this  report  are  a 
preliminary  attempt  to  understand  the  mechanisms  that 
cause  field  aligned  accelerations  in  a  H.A.N.E.  The 
1-D  and  2  1/2-D  codes  were  found  to  agree  well  with 
known  quasi linear  and  nonlinear  theory  where  checks  were 
possible.  It  was  found  in  this  reasearch  that  the  standard 
loss  cone  distribution  represented  by  a  thermalized 
unmagnetized  ion  plasma  with  a  hole  in  parallel  velocity 
space  does  not  provide  the  scattering  parallel  to  the 
magnetic  field  which  was  predicted  by  Clark  and 
Papadopoulos . *  In  the  2  1/2-D  runs  diffusion  did  fill  the 
loss  cone  but  did  not  generate  the  fluctuations  necessary 
for  diffusion  to  high  energy.  It  appears  from  comparing 
the  1-D  results  with  the  2  1/2-D  results  that  in  a 
thermalized  loss  cone  not  enough  free  energy  is  available 
to  produce  sufficiently  large  potential  fluctuations. 

This  implies  that  the  traditional  or  standard  loss  cone 
instability  produced  by  a  simple  hole  in  a  thermalized 
ion  distribution  will  not  be  energetic  enough  to  produce 
a  hard  spectrum  with  velocities  comparable  or  above  the 
expansion  velocity. 

These  results  do  not  rule  out  the  explanation  put 
forward  by  Clark  and  Papadopoulos  but  do  require  that 
their  explanation  be  more  broadly  interpreted.  As 
mentioned  in  Section  III,  there  is  a  broader  interpreta¬ 
tion  and  use  of  the  term  loss  cone.  The  term  loss  cone 
instability  has  been  applied  to  beam  plasma  type 
instabilities  where  no  particle  loss  mechanism  exists. 


{ 
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In  this  light,  there  is  still  the  possibility  that  a 

loss  cone  type  instability  might  explain  the  high  energy 

11  12 

spectrum  detected  in  the  Southern  Conjugate.  '  A 
more  appropriate  description  would  be  a  beam  plasma 
instability  of  which  there  are  several  possible.  Recall 
from  Figure  3  that  the  less  dense  plasma  distribution  was 
diffused  to  energies  ^  50%  above  its  mean  energy.  Recall  also 
that  at  15  keV  an  oxygen  atom  has  only  3  x  10  cm/sec 

g 

velocity  while  the  streaming  velocity  is  on  order  1-2  x  lo  . 
From  the  evolution  of  the  2-D  contour  plots,  it  is  con¬ 
ceivable  that  isotropization  of  the  distorted  beam  into 
parallel  space  could  produce  the  spectrum  predicted  if  one 
allows  for  higher  perpendicular  energy  at  the  outset. 

This  possiblility  will  be  examined  in  the  near  future. 

In  addition  the  inverse  mirror  geometry  of  the  magnetic 
field  will  also  be  considered  as  a  possible  acceleration 
mechanism.  Here  the  ions  must  be  magnetized. 

There  are  several  differences  between  this 
research  and  some  of  the  theory  used  to  describe  the 
loss  cone  instability.  The  major  one  occurs  because  a 
symmetry  perpendicular  to  the  magnetic  field  is  assumed 
in  most  analyses.  However,  in  the  2  1/2-D  simulations 
such  a  symmetry  does  not  exist.  The  particles  can 
gyrate  in  the  third  velocity  but  they  cannot  generate 
a  wave  spectrum  or  damp  modes  in  that  spatial  dimension. 

This  prevents  damping  of  modes  through  the  relation, 
a)  =  v^  cos  e  where  6  is  the  angle  in  the  perpendicular 
plane.  Such  deunping  could  conceivably  produce  more 
heating.  However,  if  scattering  caused  by  large  potential 
fluctuations  is  the  major  phenomena  to  be  investigated 
such  damping  would  only  inhibit  the  scattering.  An 
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example  of  such  scattering  although  caused  by  numerics  is 
seen  in  Figure  1. 

In  the  H.A.N.E.  case  such  symmetries  are  not  likely 
at  the  early  times  when  the  plasma  has  its  maximum  streaming 
energy.  As  shown  in  the  1-D  simulations,  diffusion  occurs 
on  time  scales  of  the  ion  plasma  frequency.  Using  parameters 
typical  of  Starfish  for  oxygen  and  allowing  for  a  modest 

g 

amount  of  field  compression,  .84G,  one  gets  that  w  *\»  10  , 

5  7  2  pe 

u  .  *  5  x  10  f  ai  ^  10  and  that  u>  .  ^  5  *  10  ,  The 

pX  C6  OX 

three  order  of  magnitude  difference  between  the  ion 

frequencies  implies  that  if  beaming  instabilities  do  operate 

in  the  interaction  region  they  will  occur  to  the  exclusion 

of  magnetic  behavior  on  the  part  of  the  ions.  This  result 

can  be  altered  if  the  ions  become  highly  ionized  and/or  the 

magnetic  field  is  compressed  by  amounts  approaching  an  order 

of  magnitude.  Both  events  are  possible;  in  some  Simula- 
6  18 

tions  '  magnetic  compression  of  a  factor  of  20  have  been 
produced.  Because  of  the  rapid  diffusion  one  also  has  to 
consider  the  ability  of  the  interaction  region  to  maintain  a 
loss  cone.  As  was  seen  in  the  replenishment  simulations, 
the  rate  plasma  will  stream  into  this  region  will  be  slower 
than  the  diffusion  rate.  This  implies  that  spatial  effects 
may  become  very  important  and  that  the  scattering  mechanism 
may  not  be  steady  state.  The  possibility  of  sporadic 
coupling  is  suggested  by  the  behavior  of  modes  as  seen  in 
Figures  10  and  24.  The  upshot  of  this  could  be  the 
production  of  structure  in  deposition  and  the  expanding 
debris. 

An  additional  difference  between  this  research  and 
the  standard  analytic  theory  is  that  typically  it  is  assumed 
that  the  real  frequency,  is  greater  than  the  growth 
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rate,  y*  Often  this  is  an  assumption  of  mathematical 

convenience  particularly  for  high  frequency  modes.  Our 

simulations  typically  find  y  >  to  at  least  at  the 

beginning  of  the  simulation.  As  the  simulation  continues 

this  relation  relaxes  and  y  <  u>  is  noted. 

r 

A  possible  deficiency  in  the  results  shown  was 
the  lack  of  duration  for  the  simulation.  The  simulations 
were  only  run  out  to  160/a)pe.  It  was  felt  however  that 
in  the  preliminary  phase  of  this  research  that  we  were 
looking  for  tendencies  of  the  plasma  to  be  scattered 
down  the  field  line.  If  such  tendencies  were  seen  in 
any  of  the  runs  made,  we  were  prepared  to  carry  the  run 
out.  However,  in  all  cases  the  dynamics  and  sources  of 
free  energy  had  been  diffused  by  the  instabilities  to 
levels  that  appeared  to  offer  small  chance  of  scattering. 
In  the  two  dimensional  case  none  of  the  runs  showed  any 
appreciable  electron  heating  attesting  to  the  low  level 
of  fluctuation.  In  the  parallel  direction  the  wavelength 
had  been  selected  to  miss  the  electron  Landau  damping 
region.  In  perpendicular  space  the  electron  distribution 
function  and  its  moments  changed  little. 

There  is  an  important  question  that  must  be 
addressed.  And  that  is  the  question  of  timing.  If  the 
scattering  is  occurring  as  a  direct  result  of  electro¬ 
static  instabilities  the  effect  of  the  magnetic  field 
may  be  minimal.  As  seen  in  Figures  3  and  4  only  a  small 
fraction  of  the  plasma  was  diffused  to  high  energy  but 
it  carried  most  of  the  energy  flux  in  the  center  of  mass 
frame.  The  small  density  distribution  was  heavily 
diffused  while  the  massive  one  was  not.  If  the  beams  had 
been  of  equal  mass  each  beam  would  be  diffused  but  not  to 
as  high  an  energy.  This  is  because  the  center  of  mass 


velocity  of  either  beam  relative  to  the  thermal  velocity 
would  have  been  lower.  If  the  scattering  persists  or 
occurs  for  time  greater  than  ojc^  then  the  magnetic 
field,  its  gradient  and  such  may  play  a  prominant  role 
because  a  high  beta  situation  exists.  It  then  becomes 
necessary  to  look  into  the  effects  of  Whistlers  and 
other  electromagnetic  phenomena  as  a  means  of  scattering. 
In  this  case  the  plasma  one  might  be  dealing  with  would 
be  a  thermalized  or  heavily  diffused  ion  distribution 
that  is  anisotropic,  nonmaxwellian  and  might  encompass 
most  of  the  plasma  in  the  interaction  region.  In  addition 
one  could  be  forced  to  consider  the  loss  down  the  field 
line  because  at  these  times  parallel  transport  must  be 
considered.  In  short,  it  appears  that  the  question  of 
timing  is  crucial  to  the  determination  of  whether  part  or 
most  of  the  plasma  in  the  interaction  region  is  involved 
in  the  dynamics  of  field  aligned  scattering. 

In  summary  we  feel  that  several  questions  have 
been  answered  by  this  work  and  that  the  principal 
questions  remain  to  be  answered.  Our  work  indicates  that 
the  standard  loss  cone  will  not  be  energetic  enough  to 
provide  the  speculated  high  energy  spectrum.  The  one 
dimensional  simulations  reaffirmed  the  rapid  diffusion 
that  can  be  expected  from  energetic  beam-plasma  inter¬ 
action.  Questions  remain  as  to  whether  or  not  the 
diffusion  seen  in  the  1-D  simulation  will  persist  in 
2  1/2-D  and  whether  or  not  plasma  can  be  sent  down  the 
magnetic  field  lines.  Timing  becomes  important  here 
because  it  will  determine  the  role  played  by  the  magnetic 
field.  We  suggest  that  timing  may  also  determine  the 
shape  of  the  spectrum  and  the  amount  of  energy  and  plasma 
coming  down  the  field  line. 
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A  possible  scenario  for  an  energetic  parallel 

spectrum  to  occur  is  the  interaction  of  two  relatively 

cold  beams  which  produce  the  kind  of  diffusion  upward 

in  energy  seen  in  the  one  dimensional  code.  If  this  high 

energy  plasma  is  then  diffused  parallel  by  electromagnetic 

17 

effects  such  as  ion  whistlers  then  one  could  conceivably 
end  up  with  a  very  energetic  parallel  spectrum. 

Our  future  research  is  directed  toward  answering 
these  questions.  A  mix  of  electrostatic  and  electro¬ 
magnetic  simulations  will  be  undertaken  to  consider  these 
problems  and  resolve  them. 
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ATTN:  E.  Field,  Jr 

Technology 

International  Corp 

ATTN 

W.  Boqulst 

Pennsylvania  State  University 

ATTN:  Ionospheric  Rsch  Lab 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 
SRI  International 


ATTN: 

A. 

Burns 

ATTN: 

G. 

Price 

ATTN: 

R. 

Tsunoda 

ATTN: 

J. 

Vickrey 

ATTN: 

V. 

Gonzales 

ATTN: 

J. 

Petrickes 

ATTN: 

M. 

Baron 

ATTN: 

R. 

Livingston 

ATTN: 

D. 

Neil  son 

ATTN: 

D. 

McDaniels 

ATTN: 

W. 

Chesnut 

ATTN: 

G. 

Smith 

ATTN: 

W. 

Jaye 

ATTN: 

R. 

Leada brand 

ATTN: 

C. 

Rino 

DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

TRW  Electronics  &  Defense  Sector 
ATTN:  6.  Kirchner 
ATTN:  R.  Plebuch 

Utah  State  University 
Attention  Sec  Con  Ofc  for 
ATTN:  A.  Steed 
ATTN:  D.  Burt 

ATTN:  K.  Baker,  Dir  Atmos  *  Space  Sci 
ATTN:  L.  Jensen,  Elec  Eng  Dept 

VisiDyne,  Inc 

ATTN:  W.  Reldy 
ATTN:  J.  Carpenter 
ATTN:  0.  Shepard 
ATTN:  C.  Humphrey 


Toyon  Research  Corp 
ATTN:  J.  Ise 
ATTN:  J.  Garbarlno 
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